Exceptional points in the scattering continuum 
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The manifestation of exceptional points in the scattering continuum of atomic nucleus is studied 
using the real-energy continuum shell model. It is shown that low-energy exceptional points appear 
for realistic values of coupling to the continuum and, hence, could be accessible experimentally. 
Experimental signatures are proposed which include the jump by 2t7 of the elastic scattering phase 
shift and a salient energy dependence of cross-sections in the vicinity of the exceptional point. 
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INTRODUCTION 

The structure of loosely bound and unbound nuclei is 
strongly impacted by many-body correlations and non- 
perturbative coupling to the external environment of 
scattering states and decay channels [l|, Q ■ This is partic- 
ularly important in exotic nuclei where new phenomena, 
at the borderline of nuclear structure and nuclear reac- 
tions, are expected. Some of them, like the halos [1| , the 
segregation of time scales in the context of non-Hermitian 
Hamiltonians the alignment of near-threshold states 
with decay channels [1] , and the resonance crossings 0, [3| 
appear in various open mesoscopic systems. Their uni- 
versality is the consequence of the non-Hermitian nature 
of an eigenvalue problem in open quantum systems. 

Resonances are commonly found in quantum systems, 
independently of their interactions, building blocks and 
energy scales involved. Much interest is concentrated on 
resonance degeneracies, the so-called exceptional points 
(EPs) 0. Their connection to avoided crossings and 
spectral properties of Hermitian systems [1, Q as well as 
the associated geometric phases have been discussed in 
simple models in considerable detail • The interesting 
question is their manifestation in nuclear scattering ex- 
periments. Here^ much studied case was the 2+ doublet 
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13, |13|, [IJ, |15|. Based on this example, von 
Brentano [16( discussed the width attraction for mixed 
resonances, and Hernandez and Mondragon [l3] showed 
that the true crossing of resonances can be obtained by 
the variation of two parameters in the Jordan block of 
rank two. In this latter analysis, it was shown that the 
resonating part of the scattering matrix (S-matrix) for 
one open channel and two internal states is compatible 
with the two-level formula of the R-matrix theory used 
in the experimental analysis of excitation functions of 
elastic scattering '*He(a, ao)''He [15[ and, hence, the 2+ 
doublet in ^Be may actually be close to the true reso- 
nance degeneracy. 

Properties of atomic nucleus around the continuum 
threshold change rapidly with the nucleon number, the 
excitation energy and the coupling to the environment of 



scattering states. A consistent description of the inter- 
play between scattering and resonant states requires an 
open system formulation of the nuclear shell model (see 
[E ^ U^l for recent reviews). The real-energy contin- 



uum shell model 1^ S |21|| provides a suitable unified 
framework with the help of an effective non-Hermitian 
Hamiltonian. In this work, for the first time we focus 
on a realistic model of an unbound atomic nucleus to see 
whether one or more EPs can appear in the low energy 
continuum for sensible parameters of the open quantum 
system Hamiltonian. In particular, we discuss possible 
experimental signatures of the EPs and show the evolu- 
tion of these signatures in the vicinity of the EP. Finally, 
on the example of spectroscopic factors we demonstrate 
the entanglement of resonance wave functions close to 
the EP. 



FORMULATION OF THE CONTINUUM SHELL 
MODEL 

Let us briefly review the Shell Model Embedded in the 
Continuum (SMEC) [2lj, which is a recent realization of 
the real-energy continuum shell model. The total func- 
tion space of an A— particle system consists of the set of 
square-integrable functions Q = {fpf}, used in the stan- 
dard nuclear Shell Model (SM), and the set of embedding 
scattering states V = {CI;}- These two sets are obtained 
by solving the Schrodinger equation, separately for dis- 
crete (SM) states (the closed quantum system) and for 
scattering states (the environment). Decay channels 'c' 
are determined by the motion of an unbound particle in a 
state Ij relative to the A — 1 nucleus with all nucleons on 
bounded single-particle (s.p.) orbits in the SM eigenstate 



Using these function sets, one defines projection 



operators: 
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and projected Hamiltonians: QHQ = Hqq, PHP = 
Hpp, QHP = Hqp, PHQ = HpQ. Assuming Q + V = 



2 



X, one can determine the third set of functions "''''} 
which contains the continuation of any SM eigenfunction 
tp^ in V, and then construct the complete solution in 
Q + P Recently, this approach has been extended 
to describe the two-proton radioactivity with the two- 
particle continuum [22 1. 

Open quantum system solutions in Q, which include 
couplings to the environment of scattering states and 
decay channels, are obtained by solving the eigenvalue 
problem for the energy-dependent effective Hamiltonian: 



where Hqq is the closed system Hamiltonian, c'^p ^E) 
is the Green function for the motion of a single nucleon 
in V subspace and E is the energy of this nucleon (the 
scattering energy). Index in G^p^ stands for the out- 
going boundary in the scattering problem. Tigq is non- 
Hermitian for unbound states and its eigenstates |<&q) are 
linear combinations of SM eigenstates The eigen- 

states of Hqq are biorthogonal; the left |$q,} and right 
$a) eigenstates have the wave functions related by the 
complex conjugation. The orthonormality condition in 
the biorthogonal basis reads: ($a|$/3} — Sa,p. Sim- 
ilarly, the matrix element of an operator O is Oap = 
($alO|$a). 

The scattering function 'if'^^ is a solution of a 
Schrodinger equation in the total function space: 



where 



and 



= (<P^\Hqp\Ce)/{E-£^] 



|.„ = (l + G^+'i/pQ)$„ 



Inside of an interaction region, the dominant contribu- 
tions to ^fl; are given by eigenfunctions of the effec- 
tive non-Hermitian Hamiltonian [ij: 



For bounds states, eigenvalues £a{E) of 'Hqq{E) are real 
and £a{E) = E. For unbound states, physical resonances 
can be identified with the narrow poles of the S-matrix 
0, [2^ , or using the Breit-Wigner approach which leads 
to a fixed-point condition [3, EH: 



Eo^ = Re{£^{E)) \e=e^ ; r„ = -2lm{£^{E)) \e=eA'^) 

Here it is assumed that the origin of Re {£) is fixed at 
the lowest particle emission threshold. 

An EP is a generic phenomenon in Hamiltonian sys- 
tems. In our case, the EP can appear as a result of the 



continuum-coupling term HqpG^p\e)Hpq for energies 
above the first particle emission threshold {E > 0). The 
eigenvalue degeneracies are indicated by common roots 
of the two equations Q: 



d£ 



det [Hqq {E; Vo) -£I]=0 



,y^O,l (2) 



Single-root solutions of Eq. (2) correspond to EPs as- 
sociated with decaying states. The maximal number of 
those roots is Mmax = 'n{'n — 1), where n is the number 
of states of given angular momentum J and parity tt. In 
quantum integrable models with at least two parameter- 
dependent integrals of motion one finds also double-root 
solutions which correspond to non-singular crossing of 
two levels with two different wave functions. Hence, the 
actual number of EPs in these systems is always smaller 
than Mraax ■ 

The position of EPs in the spectrum of eigenvalues of 
TYqq depends both on the chosen interaction and the 
energy E of the system. In general, eigenvalues of the 
energy-dependent effective Hamiltonian TLqq(E) need 
not satisfy the fixed-point condition ([1]) and hence need 
not correspond to poles of the S-matrix (resonances). In 
the following, we shall consider uniquely the case where 
EPs are identical with double-poles of the S-matrix. 



EXCEPTIONAL POINTS IN THE SCATTERING 
CONTINUUM OF ^"^Ne 

Let us investigate properties of EPs on the example 
of ^^Ne. SM eigenstates in this nucleus correspond to 
a complicated mixture of configurations associated with 
the dynamics of the ^^O core. Our goal is to see if EPs 
can be possibly found in the scattering continuum of 
atomic nucleus at low excitation energies and for physi- 
cal strength of the continuum coupling. SMEC calcula- 
tions are performed in ^^5/2; ■Si/2_inodel space. For 
Hqq we take the ZBM Hamiltonian [25| which correctly 
describes the configuration mixing around N = Z = 8 
shell closure. The residual coupling Hqp between Q and 
the embedding continuum V is generated by the contact 
force: Hqp = Hpq = VoS{ri — r2). For each J'^, the 
SM states \ilJi{J'^)) of the closed quantum system are 
interconnected via the coupling to common decay chan- 
nels [15f(/s:'^)op,J^; with K"" = l/2+,5/2+, and 1/2" 
which have the thresholds at E' = (the elastic chan- 
nel), 0.67 MeV, and 2.26 MeV, respectively In the ZBM 
model space, these are all possible one-proton (Ip) decay 
channels in ^^Ne. 

The size of a non-Hermitian correction to Hqq de- 
pends on two real parameters: the strength Vq of the con- 
tinuum coupling in Hqp {Hpq) and the system energy 
E. The range of relevant Vq values can be determined, 
for example, by fitting decay widths of the lowest states 
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in F. For the present Hamiltonian, experimental de- 
cay widths of the ground state 1/2^ and the first excited 
state 5/2f in ^^F are reproduced using Vq = -3500 ±450 
MeV-fm^ and Vq = -1100 ±50 MeV-fm^, respectively. 
The error bars in Vq reflect experimental uncertainties of 
those widths. The weak dependence of Ip decay widths 
on the sign of Vq is generated by the channel- channel 
coupling and disappears in a single-channel case. 
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FIG. 1: The map of J'^ = 1~ exceptional points in the con- 
tinuum of ^®Nc as found in SMEC. For more details, see the 
description in the text. 



Fig. 1 shows energies E and strengths Vq which cor- 
respond to J'^ = 1~ EPs in the scattering continuum 
of i^Ne. Decay channels [^^F^K"") ® p,J^7 with K"" = 
l/2+,5/2+, and l/2~ have been included with proton 
partial waves: Pi/2,P3/2 for K"" = 1/2+ , P3/2, /5/2, /7/2 
for K"" = 5/2+, and 51/2,^3/2 for K'' = 1/2". The 
number of 1~ SM states is 3 and, hence, the maximal 
number of 1" EPs in SMEC could be 6. Indeed, all of 
them exist at i? < 20 MeV in a physical range of Vq 
values (1100 MeV-fm^ < \Vo\ < 3500 MeV-fm^). They 
have been found by scanning the energy dependence of all 
eigenvalues over a certain range of Vq, searching for all 
real-energy crossings or width crossings (avoided cross- 
ings). Once found, we have tuned Vq to find out whether 
these crossings evolve into EPs at some combination of 
Vq and E. One should stress that the passage through 
EP always occms if, e.g., the real-energy crossing moves 
towards E = 0. Since such a crossing cannot move into 
the region E < 0, therefore it converts into an avoided 
crossing via the formation of an EP. 

The lowest EP in Fig. 1 is seen at V^"^ = -1617.4 
MeV-fm^ and E = 2.33 MeV. This EP corresponds to 
a degeneracy of the first two 1~ eigenvalues of Hqq for 
Vo < 0. Energy Ei and width Fi of lj~ and 1^ eigenvalues 
are shown in Fig. 2 as a function of the scattering en- 
ergy. For E > 2.33 MeV, width of these two eigenvalues 
grow apart very fast. Ei{E) (solid line) and E2(E) (dot- 
ted line) cross again for E ~ 3.2 MeV. At this energy, 
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FIG. 2: The upper plot exhibits the elastic scattering phase 

shifts 5piy2 (dashed-dotted line) and Sp.^^^ (dashed line) for 
p_j_i5p reaction in 1~ partial waves at around the EP (the 
double-pole of the S-matrix) with — 1^. Lower plots show 
real and imaginary parts of If (solid line) and Ij (dotted 
line) eigenvalues of the effective Hamiltonian Ti.QQ{E) as a 
function of the scattering energy E. For other details, see the 
description in the text. 



Fi and F2 are different and, hence, the corresponding 

eigenfunctions arc different as well. 

The upper part of Fig. 2 shows the phase shifts 61^ for 
p+^^F elastic scattering as a function of the proton en- 
ergy for pi/2 (dashcd-dottcd line) and P3/2 (dashed line) 
partial waves. In the partial wave Pi/2, the elastic scat- 
tering phase shift exhibits a jump by 27r at the EP with 
J'^ — 1~. This unusual jump in the clastic scattering 
phase shift is an unmistakable and robust signal of a 
double-pole of the S-matrix (EP) which persists also in 
its neighborhood, as shall be discussed below. 

Fig. 3 shows the clastic and inelastic cross sections for 
"'^^F(p, p') in the vicinity of an EP. The solid line rep- 
resents a sum of different partial contributions of both 
parities with J < 5 whereas the dashed line shows the 
resonance part of 1^ contribution in these cross sections. 
The cross sections are plotted as a function of center of 
mass scattering energy for V^""^ = —1617.4 MeV-fm^. 
The elastic cross section at the EP shows a character- 
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FIG. 3: Elastic and inelastic cross-sections in the reaction 
^^F(p, p') as a function of the proton energy E at around 
the EP (the double-pole of the S-matrix) with J'^ = \~ for 
1~ resonances only (dashed line) and for all resonances with 
J < 5 (solid line). For more details, see the description in the 
text. 



istic double-hump shape [26| with asymmetric tails in 
energy. The inelastic cross section in this case exhibits a 
single peak. Both inelastic channels [^^F(5/2+) OpiJ]j/ 

and p^F(l/2~) ®PiJ]j/ are opened at the EP. Substan- 
tial background contribution to both cross sections comes 
from broad resonances, mainly 0+ and 2+. A sharp peak 
at £■ ~ 1.65 MeV corresponds to an ordinary resonance 
2". 

The above discussion of the double-poles of the S- 
matrix (EPs) and their manifestation in the many-body 
scattering continuum concerns 1~ states. The same anal- 
ysis for J'^ = 0"'",2+ states of -'^^Ne gives qualitatively 
similar results. Also in these two cases, the number of 
EPs is maximal but only a fraction of them appears in 
the relevant range of E and Vo values. 



Behavior of scattering wave functions in the vicinity 
of the exceptional point 

A true crossing of two resonant states is accidental and, 
hence, improbable in nuclear scattering experimentation. 
In this section, we will investigate the behavior of scat- 
tering states in the vicinity of an EP (the double-pole of 
the S-matrix) as the observation of such a situation is 
more plausible. 

Fig. 4 exhibits the phase shifts Si^ for p-|-^^F elas- 
tic scattering as a function of the proton energy for 
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FIG. 4: The elastic scattering phase shifts Sp-^^^ for p-|-^^F 
reaction in 1~ partial waves at around the EP (the double- 
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-1617.4 



pole of the S-matrix) with = 1~ at Vq 
MeV-fm'^ (solid line) . Different curves correspond to different 
strength Vo of the continuum coupling: Vo=-1800 MeV-fm'^ 
(long-dashed line), -1700 MeV-fm'^ (dashed-dotted line), - 
1500 MeV-fm^ (short-dashed line) and -1430 MeV-fm^ (dotted 
line) . 



various values of the strength Vq (Vb=-1800 MeV-fm^ 
(long-dashed line), -1700 MeV-fm^ (dashed-dotted line), 
-1617.4 MeV-fm3 (solid line), -1500 MeV-fm^ (short- 
dashed hne) and -1430 MeV-fm^ (dotted line)) of the 
residual coupling Hqp = Hpq = VoS{ri — between 
Q and V subspaces. The characteristic change by a 27r 
of the elastic phase shift is seen in a broad interval -1800 
MeV-fm^ <Vq< -1500 MeV-fm^ of the continuum cou- 
pling strength. 

Fig. 5 and 6 show energies Ei and widths Fi oi 1^ 
and I2 eigenvalues as a function of the scattering energy 
for two values of Vq: -1560 MeV-fm^ (Fig. 5) and -1680 
MeV-fm3 (Fig. 6). 

The case shown in Fig. 5 corresponds to a subcritical 
coupling where two resonances cross freely in energy and 
repel in width 27|. In this regime, the scattering energy 
E corresponding to the closest approach of 1~ eigenval- 
ues in the complex plane {E ~ 2.47 MeV) is higher than 
the scattering energy corresponding to the EP at a crit- 
ical coupling V^''^'^'' =-1617.4 MeV-fm'^. Nevertheless, the 
elastic scattering phase shift shows the jump by 27r at 
the position of the EP and not at the point of the closest 
approach of eigenvalues. 

Fig. 6 shows the situation corresponding to an over- 
critical coupling where two resonances exhibit level repul- 
sion in energy and a free crossing of their widths 27| . In 
this case, the point of the closest approach of 1~ eigen- 
values in the complex plane is found at the scattering 
energy {E = 2.13 MeV) which is lower than than the 
corresponding energy for the EP. Again, the elastic scat- 
tering phase shift shows the jump by 2tt at the position 
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FIG. 5: The same as in Fig. [2] but in the subcritical regime 
of coupling (Vo = —1560 MeV-fm"'). For more details, see the 
caption of Fig. [2] and the description in the text. 
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FIG. 6: The same as in Fig. [2] but in the overcritical regime 
of coupling (Vo = —1680 MeV-fm'^). For more details, see the 
caption of Fig. [2] and the description in the text. 



of the double-pole. 

From these two examples, one can see that the charac- 
teristic jump by 27r of the elastic scattering phase shift re- 
mains a robust signature of the EP in all close-to-critical 
regimes of the coupling to the continuum: the subcriti- 
cal couphng (I Vol < iVjj'^'^^ I), the critical coupling (|Vo| = 

Vq^'^'^^I); and the overcritical coupling (|Vb| > jV'o'^'^l)' 
where real and/or imaginary parts of two eigenvalues co- 
incide. 

Next two figures show the elastic and inelastic cross 
sections for "'^^F(p, p') in the vicinity of the EP with J'^ = 
1~ in the subcritical (Fig. 7) and overcritical (Fig. 8) 
regimes of the continuum coupling. The curves shown 
by solid lines in Figs. 7,8 represent a sum of different 
partial contributions of both parities with J < 5. The 
curves shown by dashed lines exhibit the resonance part 
of 1~ contribution in these cross sections. The qualitative 
features of the cross sections for the subcritical (Vq — 
-1560MeV-fm3) and overcritical (Vb = -1680 MeV-fm^) 
couplings remain same as for the critical coupling (see 
Fig. 3). In both cases, one see a double-hump shape 
in the elastic cross sections and a single-hump shape in 
the inelastic cross section. One observes also a strong 
asymmetry in widths and heights of two peaks and a 



small shift of the position of the interference minimum in 
between the two peaks with respect to the energy which 
the EP is found for a critical coupling. 

Entangled eigenstates of the effective Hamiltonian 

Complex and biorthogonal eigenstates of the effective 
non-Hcrmitian Hamiltonian provide a convenient basis in 
which the resonant part of the scattering function can be 
expressed. These eigenstates are obtained by an orthog- 
onal and, iiigeneral, non-unitary transformation of SM 
eigenstates \y\ which is a consequence of their mixing via 
coupling to common decay channels. The same coupling 
is responsible for the entanglement of two eigenstates in- 
volved in building of an EP, as illustrated in Fig. 9 on 
the example of spectroscopic factors. 

Fig. 9 exhibits the real part of the spectroscopic fac- 
tor Re(^2)^j^g(^^i6Ne(i-)|[i5F(l/2+) p(0pi/2)]'" )') 
in ^^Ne in three regimes of continuum coupling: (a) 
the subcritical regime (Vg = —1560 MeV-fm'^), (b) the 
critical regime (vj"^ = -1617.4 MeV-fm^), and (c) the 
overcritical regime (Vb = —1680 MeV-fm'^). The solid 
(short-dashed) lines show the spectroscopic factors for 
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FIG. 7: The same as in Fig. |3] but in the subcritical regime 
of coupling (Vo = —1560 MeV-fm"^). For more details, see the 
caption of Fig. [2] and the description in the text. 
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FIG. 8: The same as in Fig. |3]but in the overcritical regime 
of couphng (Vb = —1680 MeV-fm"^). For more details, see the 
caption of Fig. [2] and the description in the text. 



^{1^){E) (<i>(l2")(i?)) eigenvalues of the effective Hamil- 
tonian Hqq{E) as a function of the scattering energy E. 
For a critical coupling (plot (b)), the spectroscopic fac- 
tors for $(1^) and $(1^) wavefunctions diverge at the 
EP (the double-pole of the S-matrix) but their sum (long- 
dashed line in Fig. 9) remains finite and constant over a 
whole region of scattering energies surrounding the EP. 
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FIG. 9: pi/2-spectroscopic factor (^''Ne(l;;)|[^''''F(l/2+) (g) 

p(0pi/2)]^ ) for ir a-nd 1^ eigenvalues of the effective Hamil- 
tonian at around the double-point of the S-matrix. For more 
details, see the discussion in the text. 



In that sense, $(!]') a-nd $(l2^) resonance wavefunctions 
form an inseparable doublet of eigenfunctions with entan- 
gled spectroscopic factors. This entanglement is a direct 
consequence of the energy dependence of coefficients baf. 

i 

in a decomposition of TYqq {E) eigenstates in the basis of 
SM eigenstates. 

One may notice that the energy dependence of Re(S'^) 
in the vicinity of the double-pole for 1^ and 1^ eigen- 
states is quite different in all three regimes of the contin- 
uum coupling. In particular, in the overcritical regime 
of coupling, an EP yields entangled states in a broad 
range of scattering energies. The strongest entanglement 
is found at the scattering energy which corresponds to 
the point of the closest approach of eigenvalues in the 
complex plane for all regimes of coupling. Obviously, the 
entanglement of resonance eigenfunctions in the vicinity 
of an EP is a generic phenomenon in open quantum sys- 
tems which is manifested in matrix elements and expec- 
tation values for any operator which does not commute 
with the Hamiltonian. 
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CONCLUSIONS 

In conclusion, we have shown in SMEC studies of the 
one-nucleon continuum that EPs exist for reahstic val- 
ues of the continuum coupling strength. In the stud- 
ied case of ^^Ne, few of those EPs appear at sufSciently 
low excitation energies to be seen in the excitation func- 
tion as individual peaks associated with a jump by 2tt 
of the elastic scattering phase shift. The occurrence of 
an EP leaves also characteristic imprints in its neighbor- 
hood, i.e. for avoided crossing of resonances. In all close- 
to-critical regimes of the continuum coupling where real 
and/or imaginary parts of the two eigenvalues coincide, 
one finds qualitatively similar features of the elastic scat- 
tering phase shift and the elastic cross-section as found 
for the critical coupling at around the EP (the double- 
pole of the S-matrix) . This gives a real chance that EPs 
or their traces may actually be searched for experimen- 
tally in the atomic nucleus. The well-known case of 2+ 
doublet in ®Be, where resonance energies and widths are 
16623±3 keV, 107±0.5 keV and 16925±3 keV, 74.4±0.4 
keV, respectively [l^ , nearly satisfies the resonance con- 
ditions in the close-to-critical regime of couplings. Var- 
ious situations in this regime have been studied experi- 
mentally in the microwave cavity 27 1. 

Avoided crossing of two resonances with the same 
quantum numbers provide the valuable information 
about the configuration mixing in open quantum sys- 
tems. As the formation of any EP in the scattering con- 
tinuum depends on a subtle interplay between internal 
Hamiltonian (Hqq) and the coupling to the external en- 
vironment of decay channels, its finding provides a strin- 
gent test of an effective nucleon-nucleon interaction and 
the configuration mixing in the open quantum system 
regime. Such tests are crucial for a quantitative descrip- 
tion of atomic nuclei in the vicinity of drip lines. 

We wish to thank W. Nazarewicz for stimulating dis- 
cussions and suggestions. 
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